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Abstract 

We show that the HERA data for the inclusive structure function F2(x,Q 2 ) for 
x < 1CP 2 and 0.045 < Q 2 < 45 GeV 2 can be well described within the color dipole 
picture, with a simple analytic expression for the dipole-proton scattering amplitude, 
which is an approximate solution to the non-linear evolution equations in QCD. For 
dipole sizes less than the inverse saturation momentum 1/Q s (x), the scattering 
amplitude is the solution to the BFKL equation in the vicinity of the saturation 
line. It exhibits geometric scaling and scaling violations by the diffusion term. For 
dipole sizes larger than 1/Q s (x), the scattering amplitude saturates to one. The 
fit involves three parameters: the proton radius R, the value xq of x at which the 
saturation scale Q s equals lGeV, and the logarithmic derivative of the saturation 
momentum A. The value of A extracted from the fit turns out to be consistent with 
a recent calculation using the next-to-leading order BFKL formalism. 



A main source of theoretical excitement about the high energy, or "small- 
x" , regime of deep inelastic lepton-hadron scattering (DIS), as currently in- 
vestigated at HERA [1], is the possibility to reach a new regime of QCD 
[2,3,4,5,6,7,8], which is characterized by high parton densities, but remains 
within the realm of perturbation theory, because the high density introduces 
a hard scale for the running of the QCD coupling a s . For sufficiently high 
energies, perturbative QCD consistently predicts that the small-x gluons in a 
hadron wavefunction should form a Color Glass Condensate (CGC) [7]. This 
is a state characterized by a hard saturation momentum Q s (x) which grows 
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rapidly with 1/x [2,5,9,10,11], and by large occupation numbers, of order l/a s , 
for all gluonic modes with transverse momentum less than or equal to Q s [5,8]. 
The saturation momentum is a measure of the gluon density in the impact 
parameter space, and thus a natural scale for evaluating a s at high energies. 

Perturbative QCD also allows us to construct non-linear evolution equations 
[6,7] (see also Refs. [2,3] for early versions of these equations) which describe 
the formation and the properties of the CGC, and its consequences for high 
energy scattering. These equations resum all powers of a s ln(l/x), together 
with the non-linear (higher- twist) effects responsible for gluon saturation and 
for the unitarization of the scattering amplitude at fixed impact parameter. 
So far, the general non-linear equations are known only to leading order (LO) 
in a s , but the accuracy is higher for the linear equation which governs the 
approach towards saturation (say, when decreasing x at fixed transverse reso- 
lution Q 2 > Ql(x)) : This is the BFKL equation [12], which is presently known 
to next-to-leading order (NLO) in a s [13], and whose renormalization-group 
improved version [14] has been recently used for a calculation of the energy 
dependence of the saturation scale [11]. The kinematical window for the valid- 
ity of the BFKL equation in the presence of saturation has been estimated as 
Q 2 (x) < Q 2 < Qs(x)/Aq CB [9]. Note that, with decreasing x, this kinematical 
region is pushed to higher Q 2 , where perturbation theory is expected to work 
better. 

It is natural to ask whether this theoretical picture can be tested against 
thef HERA data [1]. The answer is a priori not obvious: First, it is not clear 
whether the values of x at HERA are small enough for the corresponding 
saturation scale to be much larger than Aqcd- Second, the experimental points 
at HERA are correlated in such a way that the smallest values of x correspond 
also to rather small values of Q 2 (of order 1 GeV 2 , or less), where the use of 
perturbation theory becomes questionable. 

A first hint towards saturation effects at HERA came from the success of the 
simple "saturation model" by Golec-Biernat and Wusthoff (GBW) [16] which 
provided a reasonable description of the (old) HERA data [17] for x < 10~ 2 
and all Q 2 . This model has been implemented within the color dipole picture 
which is also the framework for the present analysis. The dipole picture [19] 
is a factorization scheme for DIS, which is valid at small x and is particularly 
convenient for the inclusion of unitarity corrections. Specifically, the scattering 
between the virtual photon 7* and the proton is seen as the dissociation of 
7* into a quark-antiquark pair (the "color dipole" ) followed by the interaction 
of this dipole with the color fields in the proton. This leads to the following 
expressions for the j*p cross-sections and the F 2 structure function: 

F 2 (x, Q 2 ) = (Q 2 /4n 2 a cm ) (a T + a L ), 
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a TtL (x,Q 2 )= / dzd 2 r\^ TjL (z,r,Q 2 )\ 2 a dipo i c (x,r). 



(1) 



Here, *&t,l are light-cone wavefunctions for 7* (with transverse, or longitudi- 
nal, polarization), computable within QED (see, e.g., Ref. [16] for explicit ex- 
pressions to lowest order in a em ). Furthermore, Odipoic^j r) is the cross-section 
for dipole-proton scattering (for a dipole of transverse size r), and encodes all 
the information about hadronic interactions (including unitarization effects). 
In Ref. [16], the dipole cross-section has been modelled as : 



where Q s {x) plays the role of the saturation momentum, parametrized as 
Q 2 (x) = (x /x) x GeV 2 . Saturation is visible in the fact that the dipole scat- 
tering amplitude J\f(x,r) = 1 — e~ r QsO)/ 4 approaches the unitarity bound 
Af = 1 for dipole sizes larger than a characteristic size 1/Q s (x), which de- 
creases when decreasing x. By fitting the three parameters o~o, x and A, GBW 
managed to give a rather good description of the old HERA data, for both the 
inclusive and the diffractive structure functions. More importantly — since, 
in the mean time, the simple parametrization in eq. (1) has been essentially 
ruled out [18] (see also below) by the advent of a new set of data [1], with a 
much higher accuracy — , the original success of the "saturation model" allows 
us to draw a few interesting lessons: 

i) The value of the saturation scale extracted from the fit is rather large 
(Q 2 (x) > 1 GeV 2 for x < 10~ 4 ), which suggests that a perturbative approach, 
including more of the dynamics of QCD, may work rather well in the regime 
where saturation effects become important. 

ii ) The inclusion of saturation, which tames the growth of F 2 with 1/x already 
at the hard scale Q 2 (x), may cure the small-x problem of the linear evolution 
equations [12,15] without the need to resort to soft, non-perturbative, physics. 

Hi ) The analysis of Ref. [16] suggested a remarkable regularity of the F 2 data, 
known as "geometric scaling", which has been subsequently confirmed by a 
model-independent analysis of the data [20]: The measured total cross-section 
o" 7 *p(x, Q 2 ) = ot + o~l shows approximate scaling as a function of the variable 
Q 2 /Q 2 {x) with Q 2 (x) oc l/x x and A ~ 0.3. This (approximate) scaling holds 
for small x < 10 -2 , and for all Q 2 < 450 GeV 2 . Whereas at low momenta Q 2 < 
Q 2 (x) such a scaling is indeed natural in the context of saturation, for larger 
Q 2 > Q 2 (x) it looks less natural, and it has been recently understood [9,10,11] 
as a consequence of BFKL evolution with saturation boundary conditions. The 
analysis in Refs. [9,10,11] predicts also scaling violations, in particular via the 
BFKL diffusion term, but so far these predictions have not been explicitly 
tested against the HERA phenomenology. 




(2) 
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These observations have motivated new attempts, better rooted in QCD, to- 
wards understanding the new HERA data [1] within the picture of saturation 
[18,21,22]. So far, these approaches have mostly focused upon improving the 
behaviour of the fit at large Q 2 , by including DGLAP dynamics [15]. Because 
of that, they led to rather successful global fits to the data (for x < 1CT 2 and 
all Q 2 ), but it is not clear to which extent such fits are sensitive to the details 
of the dynamics near, and at, saturation. In Refs. [18,21], the unitarization 
effects are included via the Glauber-like exponentiation of the leading-twist 
contribution, which is not fully consistent with the non-linear dynamics in 
QCD. In Ref. [22], the correct non-linear dynamics is included (to LO in a s ), 
together with DGLAP dynamics, by a combined numerical analysis of the 
Balitsky-Kovchegov (BK) equation [6] and of DGLAP equation. But also in 
that case, the success of the global fit turns out to be very much dependent 
upon the inclusion of DGLAP: Without the latter, that is, by using the BK 
equation alone, the fit of Ref. [22] works satisfactory only for rather low Q 2 
(up to a few GeV 2 ). Besides, being fully numerical, the analysis in Ref. [22] 
does not allow to simply test some qualitative features of BFKL dynamics 
towards saturation, like geometric scaling and its violations. 

Our aim in this Letter is to present a new analysis of the HERA data, restricted 
to the kinematical range where one expects important high density effects - 
namely, x < 10~ 2 and Q 2 < 50 GeV 2 — , with the purpose of showing that 
these data are consistent with our present understanding of BFKL evolution 
and saturation. The upper limit on Q 2 has been chosen large enough to include 
a significant number of "perturbative" data points, but low enough to justify 
the emphasis on BFKL, rather than DGLAP, evolution (for the small-x values 
of interest). In fact, 50 GeV 2 is roughly consistent with the estimated upper 
bound Q^(x) / Aqqy) of the kinematical window for BFKL behaviour [9], given 
our subsequent findings for Q s (x). Within this kinematical range, we shall 
provide a reasonable fit (x 2 /d.o.f. ~ 0.8 — 0.9) to the new HERA data for F 2 
based on a simple, analytic, formula for the dipole scattering amplitude, which 
is an approximate solution to the non-linear evolution equations in QCD [6,7]. 
We shall refer to this fit as the "CGC fit". 

It turns out that our fit involves the same set of free parameters (that is, a , 
x and A) as the GBW "saturation model". The need for the first two pa- 
rameters, o"o and xq, reflects the fact that, even in a first principle calculation 
including saturation, some aspects of the calculation remain non-perturbative: 
the impact parameter dependence of the scattering amplitude and the initial 
condition at low energy. Concerning the first aspect, we shall treat the proton 
as a homogeneous disk of radius R. (A more realistic impact parameter de- 
pendence, which requires a model, can be implemented along the lines of Refs. 
[21,22]. This is left for a later publication [23].) Then, o"di po ie(^, r) = o- Af(x, r) 
with uq = 2ixR 2 and Af(x, r) given by the solution to the homogeneous version 
of the non-linear evolution equation. Within the approximation that we shall 
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use to solve this equation, the initial condition is fully characterized by one 
parameter: the value Xq of x at which Q s equals lGeV. The third parameter, A, 
which controls the energy dependence of the saturation scale, is conceptually 
different, since this can be computed in perturbation theory [2,9,10,11]. As 
already mentioned, this parameter is presently known to NLO accuracy [11]. 
More precisely, the calculation in Ref. [11] shows that \{Y) = d\nQ 2 (Y)/dY 
(with Y = \nl/x) is not simply a constant, but rather a slowly varying func- 
tion, which decreases from A ~ 0.30 for Y = 5 to A ~ 0.27 for Y = 15, with a 
theoretical uncertainty of about 15%. Unfortunately, this uncertainty is still 
too large to permit a good description of the data: A is the exponent which 
controls the growth of F2 with Y, so the data are very sensitive to its precise 
value. To cope with that, we shall treat A as a free parameter, that we shall 
fit from the data. Remarkably, the value of A that will come out from the fit, 
namely A 0.25 — 0.29, is consistent with the theoretical prediction of Ref. 
[11], within the theoretical uncertainty alluded to above 1 . 

The function Af(Y, r) will be constructed by smoothly interpolating between 
two limiting behaviours which are analytically under control: the solution 
to the BFKL equation for small dipole sizes, r <C 1/Q s (x), and the Levin- 
Tuchin law [24] for larger dipoles, r 3> l/Q s (x). As we shall check later, the 
quality of the fits is not very sensitive to the details of the interpolation, nor 
to the precise form of the approach towards the limit J\f — 1. On the other 
hand, the data are quite sensitive to the details of the scattering amplitude 
at smaller sizes r < \/Q s (x), and thus provide a test of BFKL dynamics. 
In this range, J\f(Y, r) will be obtained via the saddle point approximation 
to the LO BFKL equation, followed by an expansion to second order around 
the saturation saddle point. As known from Refs. [9,10], the first term in 
this expansion exhibits geometric scaling, while the second, "diffusion", term 
brings in scaling violations, which turn out to be crucial for understanding 
the HERA data. The LO formalism is chosen for its simplicity: it provides an 
explicit expression for J\f(Y,r), whose physical interpretation is transparent. 
This is further motivated by the observation [11] that, with increasing Y, the 
predictions of the (RG-improved) NLO BFKL formalism get closer and closer 
to those of the LO equation with running coupling a s (fi 2 ), in which the scale 
/i 2 for running is set either by the size of the dipole (/i 2 = C/r 2 ) [10], or 
simply by the saturation momentum (/z 2 = Q 2 (x)) [9]. To account for the 
lack of accuracy of the LO formalism at non-asymptotic Y, we shall treat, as 
announced, the saturation exponent A — which is the critical parameter for 



On the other hand, this value for A is significantly smaller than that predicted 
by the LO BFKL equation with running coupling [2,9,10,11], which varies between 
0.43 and 0.37 when increasing Y from 5 to 10 [11]. This discrepancy may explain 
why the pure BK part of the fit in Ref. [22] (based on a numerical solution to the 
BK equation [6] with running coupling) is unable to accurately describe the data 
with Q 2 > 1 GeV 2 and small x. 
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describing the data — as a free parameter. 



Let us now briefly introduce our formulae. For both fixed coupling a s , or 
running coupling a s (Q s (Y)), the solution to the BFKL equation can be written 
in Mellin form as (for more details, see [9,10]): 

AT(r,r) = /^(r 2 ^) 7 e^W ( 7 ) (3) 

c 

where Q is a reference scale of order A QCD , x(t) = 2-^(1) — V'(t) — ~~ 7) 
with ^(7) = dlnT('~f)/d'j, ^0(7) is the initial condition, and the function 
h(Y) depends upon our assumption on the running of a s : For fixed coupling, 
h(Y) = a s Y with a s = a s N c /n, while for a running coupling a s (Q s (Y)), h(Y) 
is determined by dh/dY = a s (Q s (Y)), with h(0) = 0. 

In the saddle point approximation, valid when h(Y) and p = \n(l /r 2 Ql) are 
large, 

J\f(Y, r) ~ e HY)F(io(K),ii) (4) 



where TZ = p/h(Y), F{^,TZ) = — 7"7£ + x(7), and 7o(7?.) is the saddle point for 
a given 1Z, determined by: 



= -K+ x'(To(K)) = 0. (5) 

7=70 (R-) 



Note that, in the approximation above, we have neglected the initial condition 
A/o(7), as well as the effect of the Gaussian fluctuations around the saddle 
point: Indeed, these are slowly varying functions, which contribute at most 
logarithmic terms (like In ft, or lnp) to the exponent in Eq. (4). 

The saturation condition is written as: Af(Y, r) = Af for r = 1/Q S (Y), where 
Ao is a number of order one (its precise value is a matter of convention). To 
the accuracy of interest, this condition implies 

F( l0 (n),n)\ n=iis = - lo (K s )K s + X (lo(K s )) = 0, (6) 
where TZ S = Ps (Y)/h(Y) and p s (Y) = \n(Q 2 s (Y)/Q 2 ). 

Eq. (6) shows that 1Z S is a pure number. Together with Eq. (5), it allows us 
to compute both 1Z S and 7 S = 7o(7£ s ) (the saddle point in the vicinity of 
the saturation line Q 2 = Q 2 (Y)). One obtains [2,9,10] : 7 S = xilsj/x'ds) = 
0.627... and 1Z S = x'ils) — 4.883.... Note that 7 S is not the same as the usual 
BFKL saddle point 70 = 1/2 (the latter would be obtained by letting 1Z — > 
in Eq. (5)). Whereas 7 o describes the evolution with Y at fixed, and relatively 
low, Q 2 (here, Q 2 ~ 1/r 2 ), 7s corresponds rather to an evolution where, when 
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increasing Y, Q 2 is correspondingly increased, in such a way that the condition 
Q 2 ~ Q 2 S (X) ( or P ~ Ps(Y)) remains satisfied. 

The dipole sizes r that we are interested in are smaller than, but relatively 
close (in logarithmic units) to, 1/Q S (Y). It is therefore a very good approxi- 
mation to compute MiY, r) by expanding the exponent in Eq. (4) to second 
order in p - p s = \n(l /r 2 Q 2 (Y)). One then obtains [9,10]: 



where j3 = x"(ls) — 48.518.... Eq. (7) shows that, in the vicinity of the satura- 
tion line (p ~ p s ), where the quadratic term in the exponent can be neglected, 
the scattering amplitude scales as a function of r 2 Q 2 {Y). This is geometric 
scaling. The power 7 S (or, more precisely, the difference 1 — 7 S 0.37) is often 
referred to as an "anomalous dimension" . The quadratic, "diffusion" , term in 
Eq. (7) is responsible for scaling violations. 

Now that we have found 1Z S , one could use its original definition, namely 1Z S = 
(l/h(Y))\n(Q 2 (Y)/Ql), together with the equation dh/dY = a s (Q s (Y)) to 
compute Q S (Y) [9,10]. However, if one does so, then the resulting expression 
for QsiY) will not be accurate enough to describe the data. Rather, we shall 
rely on the NLO calculation in Ref. [11] to conjecture that p s (Y) = XY with A 
a pure number, to be fitted from the data. As for the other coefficients which 
appear in Eq. (7), namely 7 S and the ratio k = (3/lZ s , these will be kept as in 
the LO BFKL approximation since (a) their numerical values do not change 
appreciably when going to NLO [11], and (b) the fit is not very sensitive to 
their precise values, as we shall check. This last feature — the stability of 
the fit with respect to small changes in the values of the coefficients 7 S and k 
- is important since, given the present status of the theory, it seems to be 
quite difficult to improve our theoretical control of these parameters (e.g., by 
going to NLO accuracy). The solution J\f(Y,r) to the (RG-improved) NLO 
BFKL equation with saturation boundary condition is known [11] only for 
asymptotically large Y, and thus cannot be applied to the kinematical range 
at HERA without introducing additional free parameters. Therefore, using 
that (rather complicated) solution in a fit would obscure the physics without 
automatically improving the predictive power of the theory. 



To summarize, the dipole cross-section that we shall use in the CGC fit reads 
0di P oie(:r, r) = 2nR 2 M(rQ s , Y), with 




(7) 




M(rQ 8 ,Y) = l-e 



-a In 2 



(brQ a ) 



for rQ s > 2, 



(8) 
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where Y = ln(l/x), Q s = Q s (x) = (x /x) x / 2 GeV, and we have redefined Q s 
in such a way that J\f(rQ s , Y) = Ao for rQ s = 2 (to facilitate the comparison 
with the GBW results based on Eq. (2)). The expression in the second line 
has the correct functional form for r 3> 2/Q s , as obtained either by solving 
the BK equation [24], or from the theory of the CGC [8] (see also [25] for a 
more careful discussion). This is strictly valid only to LO accuracy, but here 
it is used merely as a convenient interpolation towards the 'black disk' limit 
A" — 1. (The details of this interpolation are unimportant for the calculation 
of <7 7 * p .) The coefficients a and b are determined uniquely from the condition 
that N(rQ s , Y) and its slope be continuous at rQ s = 2. The overall factor 
Ao in the first line of Eq. (8) is ambiguous, reflecting an ambiguity in the 
definition of Q s . We shall repeat our fit for various values of Ao between 0.5 
and 0.9 and check that the results of the fit change only little. 

As already mentioned, the coefficients 7 S and k are fixed to their LO BFKL 
values: 7 S = 0.63 and k = 9.9. We shall use the same photon wavefunctions 
\J/t,l as in Refs. [16,18,21]. These involve a sum over the active quark flavors, 
and in our fits we shall work with three quarks of equal mass m q , for which 
we shall choose three different values: m q = 140 MeV, 50 MeV, and 10 MeV. 
(Further analysis of the m^-dependence of the fit, as well as the inclusion of 
the heavier, charm quark are left for Ref. [23].) Thus, the only free parameters 
of the fit are R, x and A, as announced. 

Our fit has been performed for the F 2 data at ZEUS (the first two references 
in [1]) with x < 10~ 2 and Q 2 between 0.045 and 45 GeV 2 (156 data points). 
Only ZEUS data have been considered because there is a mismatch between 
HI and ZEUS concerning the data normalization, and it is only ZEUS which 
has data at very low Q 2 , in the saturation region. (Adding the HI data in the 
fit would have implied a normalization adjustment, and would have put more 
weight to the moderate and large Q 2 data.) Statistical and systematic errors 
have been added in quadrature. 

The values obtained for the various parameters, and the x 2 /d.o.f. for the fits 
corresponding to different choices for Ao and m q , are shown in Tables 1 and 
2. Note, in particular, the good stability of the result for A, which changes 
only by 15% (within the range A = 0.25 — 0.29) when varying Ao and m q . 
As anticipated, this value of A is in agreement with the theoretical calculation 
at NLO in Ref. [11]. In Table 1, we also show, for comparison, the results 
obtained when performing a fit to the same set of data with the GBW cross- 
section (2). Clearly, the corresponding x 2 /d.o.f. = 1.59 is relatively poor. Note 
also that the value of Q 2 s (x) emerging from the CGC fit is smaller — roughly, 
by a factor of 2 for Ao = 0.7 — than the corresponding value for the GBW 
fit (since x is correspondingly smaller). 

In Figs. 1 and 2, the results of the fit are plotted against the data for Ao = 0.7 



8 



and m q = 140 MeV. In Fig. 2, we have also shown the extrapolation of our 
fit towards larger values of x and Q 2 (outside the range of the fit), together 
with the corresponding HERA data, in order to emphasize that the deviations 
from BFKL dynamics — due notably to the presence of the valence quarks 
and to the DGLAP evolution [15] — do eventually show up, as expected. 

In both Figs. 1 and 2, we also show (with dashed line) the prediction of the 
BFKL calculation without saturation, as obtained by extrapolating the for- 
mula 2 in the first line of Eq. (8) to arbitrarily large rQ s (with the same values 
for the parameters as in Table 1, and the infrared diffusion term switched off 
for rQ s > 2). This pure BFKL fit shows a too strong increase with \jx at small 
Q 2 , as expected from similar analyses in the literature [26]. On the other hand, 
the complete fit, including saturation, works remarkably well down to the low- 
est values of Q 2 that we have included. Since, a priori, such low values of Q 2 
seem to be completely out of the reach of perturbation theory, it is important 
to explain why, in the context of saturation, it is still meaningful to approach 
these data via the previous calculation: 



A/"o / model 


0.5 


0.6 


0.7 


0.8 


0.9 


GBW 


x 2 


146.43 


129.88 


123.63 


125.61 


133.73 


243.87 


xVd.o.f 


0.96 


0.85 


0.81 


0.82 


0.87 


1.59 


x (xlO- 4 ) 


0.669 


0.435 


0.267 


0.171 


0.108 


4.45 


A 


0.252 


0.254 


0.253 


0.252 


0.250 


0.286 


R (fm) 


0.692 


0.660 


0.641 


0.627 


0.618 


0.585 



Table 1 

The CGC fits for different values of Mo and 3 quark flavors with mass m q = 140 
MeV. Also shown is the fit obtained by using the GBW model, Eq. (2). 

The reason is that, for the low Q 2 data at HERA, the associated values of x are 
also very small, so that the corresponding saturation scale is relatively large: 
Q 2 (x) ~ 1.3 - 2.3 GeV 2 for x ~ 10~ 5 - 1(T 6 . Thus, for these data, we have 
Q 2 <C Q 2 (x), which implies that the large, non-perturbative dipoles (r ^ 1/Q) 
are deeply in the saturation regime, so for them A/" ~ 1 independently of 
the detailed mechanism leading to such a strong absorbtion. Clearly, the qq 
pairs with transverse sizes r > 2 GeV -1 have a strong overlap with hadronic 
states, like the p meson, whose correct treatment goes beyond the scope of 

2 Note that this is not the standard BFKL formula, since the saddle point 7 S is 
different from 70 = 1/2. But the behaviour illustrated by the dashed line fit in Figs. 
1 and 2 is representative for the BFKL fits without saturation [26]. For instance, 
a similar behaviour can be found by extrapolating the fit in the last paper of Ref. 
[26] towards low Q 2 . 
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with N n =0.7 and m =140 Me\ 



Fig. 1. The F2 structure function in bins of Q 2 for small (upper part) and moderate 
(lower part) values of Q 2 . The experimental points are the latest published data 
from the HI and ZEUS collaborations [1]. (The HI data have been rescaled by a 
factor 1.05 which is within the normalization uncertainty.) The few data points at 
the lowest available Q 2 (0.045, 0.065 and 0.085 GeV 2 ) are not displayed although 
they are included in the fit. The full line shows the result of the CGC fit with 
Mo = 0.7 to the ZEUS data for x < 10~ 2 and Q 2 < 45 GeV 2 . The dashed line shows 
the predictions of the pure BFKL part of the fit (no saturation). 

10 




Fig. 2. The same as in Fig. 1, but for large Q 2 . Note that in the bins with 
Q 2 > 60GeV 2 , the CGC fit is extrapolated outside the range of the fit 
(Q 2 < 50 GeV 2 and x < 10 -2 ), to better emphasize its limitations. 
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m g = 50 MeV 


m q = 10 MeV 


Kfn 


U.o 


n 7 
U. ( 


n o 

u.y 
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n 7 
U. i 
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X 




lUO.tJZ 


1 ns 7fi 

1UO. / u 


1 AQ 97 


1 f!7 R/l 




y 2 /d.o.f 


0.97 


0.71 


0.71 


0.98 


0.70 


0.70 


x (xlO" 4 ) 


2.77 


0.898 


0.333 


3.32 


1.06 


0.382 


A 


0.290 


0.281 


0.274 


0.295 


0.285 


0.276 


i? (fm) 


0.604 


0.574 


0.561 


0.593 


0.566 


0.554 



Table 2 

The CGC fits for three values of A/"o and quark masses m q = 50 MeV (left) and 
m q = 10 MeV (right). 

perturbation theory. But in the calculation of the DIS cross-section a 7 » p at 
small x, all such states are summed up with the same weight (since they 
all have an amplitude N ~ 1), so this sum simply accounts for the total 
probability for hadronic components with Q 2 < 1 GeV 2 in the virtual photon 
wavefunction. It is reasonable to assume that this total probability is correctly 
given by perturbation theory (i.e., by summing over all qq states with sizes 
r > 2 GeV" 1 ), although the way how this probability is distributed among 
various meson states is truly non-perturbative 3 . Note also that, because of 
saturation, the distribution of the dipole sizes in the integrand of Eq. (1) 
when Q 2 ^ Q 2 s {x) is pushed towards smaller dipoles, at least to logarithmic 
accuracy [16,21,23] : the typical dipole sizes which contribute to the integration 
are logarithmically distributed in the range Q 2 1/r 2 Q 2 {x). 

As anticipated, we have checked that neither the quality of the fit, nor the 
value of the parameters, change appreciably when the coefficients 7 S and k in 
Eq. (8) are modified by about 10 % (but the fit appears to be more sensitive 4 
to the value of 7 S then to that of k) [23]. Such variances are representative 
for the theoretical uncertainties of our BFKL calculation, so it is reassuring 
to see that the success of our fit is not dependent upon a fine-tunning of these 
parameters. 

Finally, in order to study the role of scaling violations, we have performed a 
fit with a pure scaling function [27], namely: M(rQ s ) = Af (rQ s /2) 2 ' y (with 
7 a free parameter) for rQ s < 2, and saturation (as in the second line of 
Eq. (8)) for rQ s > 2 (see results in Table 3). We have found that, although 



3 We would like to thank Al Mueller for pointing out this argument to us. 

4 One should also keep in mind that these two parameters are actually correlated, 
so it makes no sense to consider strong, independent, variations in their values. 
For instance, in the LO BFKL formalism that we are using here, they are both 
determined by properties of the function x{l) near the saturation saddle point 7 S . 
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this fit involves 4 parameters, its quality is relatively poor: x 2 /d.o.f. ~ 1.4. 
More interestingly, we have found that the "best" value of 7 for describing the 
data in the kinematical range of interest (and within this scaling Ansatz) is 
7 ps 0.84, which is significantly smaller than one (at variance with the GBW 
"saturation model", and also with DGLAP evolution, which would predict 
a power-law behaviour with 7 = 1 up to logarithmic corrections), but also 
significantly larger than the LO BFKL value 7 S = 0.63. This clearly shows 
that the presence of the diffusion term in Eq. (7) has been essential for the 
success of the 3-parameter CGC fit: This term violates scaling for r < 2/Q s , 
and enhances the effective "anomalous dimension" from j s to 



= dlnM(rQ„Y) = K2/rQ.) 
WV" r >- dlii(4/r 2 QJ) " + k.\Y ■ [> 

The difference 7eff(rQ s ,^) — 7s decreases with Y, but increases with the de- 
viation p — p s = 2ln(2/rQ s ) from the saturation line. This behaviour is illus- 
trated in Fig. 3, where we have also included, for comparison, the pure scaling 
functions with 7 S = 0.63 and 7 = 0.84 (the average "anomalous dimension" 
preferred by the fit in Table 3). Note that, for any x, there is a significant 
range, roughly at 1 < rQ s < 2, within which the geometric scaling approxi- 
mation with 7 S = 0.63 works quite well. On the other hand, for rQ s <C 1, 7 e g- 
gets closer to one 5 , and thus mimics DGLAP behaviour; within the range of 
the fit, this is important only for the data with relatively large Q 2 and not so 
small values of x, which are more sensitive to such small-size dipoles. 



A/o 


x 2 


X 2 /d-0.f 


x (xl0~ 4 ) 


A 


R (fm) 


7 


0.7 


215.70 


1.42 


3.79 


0.313 


0.572 


0.845 



Table 3 

A 4 parameter fit for geometric scaling + saturation. 

Further variations and checks of our fits, together with more applications to 
the HERA phenomenology, will be presented somewhere else [23]. 
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5 Taken literally, our formula implies 7 e g — > 00 as r — > 0, which together with 
Eq. (8) implies that, when r — > 0, the dipole amplitude vanishes faster than any 
power of r. However, this (unphysical) behaviour has no influence on the fit, since in 
the kinematical range of interest the integrated result for almost insensitive 

to very small dipoles. 
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Fig. 3. The dipole amplitude for two values of x, compared to the pure scaling 
functions with "anomalous dimension" 7 = 7 S = 0.63 and 7 = 0.84. 
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